Preprint typeset using WT^^ style cmulatcapj v. 04/21/05 



SPITZER MID-INFRARED SPECTROSCOPY OF DISTANT X-RAY LUMINOUS AGN 

Kate Brand\ Dan W. Weedman^, Vandana Desai^, Emeric Le Floc'h**, Lee Armus'"', Arjun Dey*', Jim R. 
HouCK^, BuELL T. Jannuzi^, Howard A. Smith'', B. T. Soifer^'^ 

ABSTRACT 

We present mid-infrared spectroscopy of a sample of 16 optically faint infrared luminous galaxies 
obtained with the Infrared Spectrograph (IRS) on the Spitzer Space Telescope. These sources were 
jointly selected from Spitzer and Chandra imaging surveys in the NDWFS Bootes field and were 
selected from their bright X-ray fluxes to host luminous AGN. None of the spectra show significant 
emission from polycyclic aromatic hydrocarbons (PAHs; 6.2/im equivalent widths <0.2/im), consistent 
with their infrared emission being dominated by AGN. Nine of the X-ray sources show 9.7/im silicate 
absorption features. Their redshifts are in the range 0.9 < z < 2.6, implying infrared luminosities of 
log(L/fl) = 12.5— 13.6 Lq. The average silicate absorption strength is not as strong as that of previously 
targeted optically faint infrared luminous galaxies with similar mid-infrared luminosities implying that 
the X-ray selection favors sources behind a smaller column of Si-rich dust than non- X-ray selection. 
Seven of the X-ray sources have featureless power-law mid-IR spectra. We argue that the featureless 
spectra likely result from the sources having weak or absent silicate and PAH features rather than 
the sources lying at higher redshifts where these features are shifted out of the IRS spectral window. 
We investigate whether there are any correlations between X-ray and infrared properties and find 
that sources with silicate absorption features tend to have fainter X-ray fluxes and harder X-ray 
spectra, indicating a weak relation between the amount of silicate absorption and column density of 
X-ray-absorbing gas. 

Subject headings: galaxies: active — galaxies: starburst — infrared: galaxies — quasars: general 



1. INTRODUCTION 

Mid-infrared spectroscopy is a useful tool in investigat- 
ing the extent to which starbursts and AGN contribute 
to the luminosity of luminous infrared galaxies (LIRGs; 
Sanders et al. 1988). The infrared spectra of local LIRGs 
have been studied extensively with ISOPHOT, ISOCAM 
and SWS on the Infrared Space Telescope (ISO; e.g., 
Lutz et al. 1998; Genzel et al. 1998; Rigopoulou et al. 
1999; Tran et al. 2001) and more recently with the in- 
frared spectrograph (IRS; Houck et al. 2004) on the 
Spitzer Space Telescope (Weedman et al. 2005; Brandl 
et al. 2006; Armus et al. 2007; Desai et al. 2007). Spitzer 
IRS has also extended these studies to fainter sources at 
higher redshifts. A particularly interesting population 
of high redshift optically faint infrared luminous sources 
(Yan et al. 2004; Brand et al. 2007; Fiore et al. 2008; 
Dey et al. 2008) has been discovered using the Multiband 
Imaging Photometer (MIPS; Rieke et al. 2004). Various 
observing programs with Spitzer IRS have found that 
MIPS sources at flux density levels of fi/(24/zm) > 0.8 
mjy with optical magnitudes i? > 24 are typically at z ~ 
2-3 (Houck et al. 2005; Weedman et al. 2006b Yan et al. 
2007; Sajina et al. 2007). The infrared spectra exhibit 
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a wide range of properties from deep 9.7/im silicate ab- 
sorption features and/or strong polycyclic aromatic hy- 
drocarbon (PAH) emission features to featureless power- 
law continua. This obscured population is important in 
determining how much light originates from starbursts 
and AGN at epochs in the Universe when their luminos- 
ity density was at its maximum. Given that their optical 
faintness makes optical follow-up difficult, understand- 
ing how to characterize these sources using their infrared 
spectra is important in understanding their energetics. 

One approach in understanding the contribution of 
starbursts and AGN to the infrared emission of opti- 
cally faint infrared sources is to investigate the infrared 
spectra of sources also selected in other wavebands. For 
example, at one extreme we can study sources selected 
on the basis of large far-infrared flux densities, which we 
might expect to contain cool dust and show signatures of 
starbursts. In Brand et al. (2008), we presented the IRS 
spectra of 11 70/im-selected optically faint infrared lumi- 
nous galaxies and found that seven had large 6.2/im PAH 
equivalent widths, consistent with them being powerful 
starburst galaxies. Four of the galaxies showed deep sil- 
icate absorption features and their properties suggested 
that AGN are likely to be the dominant origin of their 
large mid-infrared luminosities. 

At the other extreme, we can consider sources selected 
on the basis of X-ray flux, which we might expect to show 
spectral signatures of AGN. Weedman et al. (2006a) pre- 
sented IRS spectra of 9 X-ray-loud optically faint infrared 
luminous galaxies selected from the Spitzer Wide-Area 
Infrared Extragalactic Survey (SWIRE; Lonsdale et al. 
2004). Eight sources have sihcate absorption features 
and one source has a featureless power-law mid-IR spec- 
trum. 

In this paper, we present Spitzer IRS spectra of 16 
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optically faint infrared galaxies in the Bootes field of the 
NOAO'Decp Wide-Field Survey (NDWFS; Jannuzi & 
Day 1999) which are selected to be X-ray loud. This sig- 
nificantly increases the number of X-ray loud optically 
faint infrared luminous galaxies thus far studied, and ex- 
tends the sample to higher X-ray and 24/im luminosities. 

A cosmology of Hq = 70 km Mpc"""^, f2M=0.3, and 
riA=0.7 is assumed throughout. 

2. SOURCE SELECTION 

Our sample is selected from the Bootes field of the 
NOAO Deep Wide-Field Survey (NDWFS; Jannuzi & 
Dey 1999). The combination of the large survey area 
(8.2 deg^ with overlapping MIPS and optical coverage) 
and existing Spitzer MIPS and IRAC, Chandra, and 
deep optical imaging data makes this an ideal data-set 
for multi-wavelength comparisons of sources. The MIPS 
data were obtained with an effective integration time at 
24/im of ^90 s per sky pixel, reaching a 5 cr detection 
limit of ^ 0.3 mJy for unresolved sources. The resulting 
source lists contain w28,000 sources with fj,(24;um)>0.3 
mJy. The XBootes survey was undertaken using ACIS-I 
on the ChandraX — rayObservatory with 5 ks extosures 
over the entire Bootes field (Murray et al. 2005; Kenter 
ct al. 2005; Brand et al. 2006a). It comprises 3293 point 
sources with 4 or more X-ray counts, corresponding to a 
limiting flux of f(o.5_7fcey)=8.1xlO~^^ ergscm"^ s . An 
8.5 deg^ region of the Bootes field has been mapped in all 
four Spitzer IRAC bands by the IRAC Shallow Survey 
(Eisenhardt et al. 2004), reaching 5a Umits of 6.4, 8.8, 
51, and 50 /xJy at 3.6, 4.5, 5.8, and 8 /zm, respectively. 

Using the XBootes and MIPS 24/xm catalogs, we se- 
lected a sample of 16 X-ray loud sources in the Bootes 
field for Spitzer IRS observations. We selected all 13 
of the sources with >4 X-ray counts (fo.5-7fcey > 8.1 x 
10-15 ergs cm-2 s'^), f^(24^m)> 2 mJy, and R > 22 
mag (the "main" sample). We supplemented this sample 
with 3 extra sources (BootesX13, 14, and 15) which were 
slightly fainter at 24 /im (fj^ (24/im) >1 mJy) but which 
were particularly powerful X-ray sources (hereafter de- 
noted as the "X-ray bright" sub-sample with 48, 34, and 
19 X-ray counts corresponding to {o.5-7keV= 9.2, 6.5, and 
3.7 X lO^i^crgs cm^^ s^^ respectively). Because of their 
different selection criteria, we identify the X-ray bright 
sub-sample separately throughout the analysis and con- 
sider the conclusions with and without their inclusion. 

The basic optical and infrared properties of the sample 
are shown in Table 1. 10/16 of the sources have fairly red 
optical to infrared colors and qualify as Dust Obscured 
Galaxies (DOGs) based in the criterion (i?-[24]>14) pre- 
sented in Dey et al. 2008. None of our sample have red- 
shifts from optical spectroscopy. The X-ray properties 
are presented in Table 2. At z > 1, the X-ray luminosi- 
ties are far larger than can be attributed to a starburst 
and must be attributed to an AGN. At z = 0.9 (the low- 
est measured source rcdshift; see below), the minimum 
X-ray flux, fo.5-7fcey = 8.1 x 10-^^ ergs cm"^ s"^ cor- 
responds to an X-ray luminosity of » 3 x 10^^ ergs s~^. 
If this was all due to star-formation, it would imply an 
unfeasibly large star- formation rate of « 5000 M© yr~-^ 
(Grimm et al. 2003). 

3. OBSERVATIONS AND DATA ANALYSIS 



The spectroscopic observations were made with the 
IRS Short Low module in order 1 only (SLl) and with 
the Long Low module in orders 1 and 2 (LLl and LL2), 
described in Houck et al. (2004). These orders give low 
resolution spectral coverage from ~8 /um to ~35 /xm. The 
integration times for individual sources are given in Ta- 
ble 3. Data were processed with version 13.0 of the SSC 
pipeline, and extraction of the source spectra was done 
with the SMART analysis package (Higdon et al. 2004). 
For details of the observational and data analysis tech- 
niques, we refer the reader to Brand et al. (2008). 

4. RESULTS 

The IRS spectra of all 16 X-ray selected sources are 
shown in Figure 1. Although the IRS spectrum of 
Houckl3 was presented in Houck et al. (2005), it sat- 
isfies the selection criteria of the main sample and we 
include it again here for completeness. 9/13 of the main 
sample (and 9/16 of the total "main -I- X-ray bright" 
sample) exhibit identifiable 9.7/im silicate absorption fea- 
tures (Sio <-0.4^; hereafter described as "absorption- 
dominated" sources and denoted by "abs" in Table 3). 
A further 4/13 or the main sample (and 7/16 of the total 
"main + X-ray bright" sample) exhibit power-law spec- 
tra with no convincing features (hereafter described as 
"power-law" sources and denoted by "pow" in Table 3); 
either they have only weak silicate absorption features 
or they are at z > 2.6 and the absorption feature does 
not fall within the observable IRS bandwidth. BootesX6 
may have a weak silicate emission feature if it is at z=0.9. 
Given that this is highly uncertain without other con- 
vincing features in the spectrum, we classify BootesX6 as 
a power-law source. None of the sources show convincing 
PAH emission features, although, given their uncertain 
or unknown redshifts, there may be weak PAH emission 
in some individual sources. We discuss each class in turn. 

4.1. Absorption-Dominated Sources 

For the absorption-dominated sources, redshifts are de- 
termined primarily from the localized maximum in the 
IRS continuum, or "hump" , which is blueward of the ab- 
sorption feature. This feature is produced by absorption 
on either side of the feature, and with a possible contri- 
bution of 7.7/im PAH emission in composite sources. We 
determine redshifts by assuming this hump has a rest- 
frame wavelength of 7.9/im (see average spectra in Hao 
et al. (2007) and Spoon et al. (2007)). Comparison of 
IRS redshifts with optical redshifts for similar absorbed 
sources (e.g.. Brand et al. 2007; Dey et al. 2008) suggests 
that they can be uncertain to ±0.2 in z. 

For the absorption-dominated sources, the overall na- 
ture of the infrared spectra can be best considered by 
assembling the average IRS spectrum. This allows con- 
sistent features within individual spectra to be seen with 
improved signal-to-noise. This was achieved by wave- 
length correcting each individual unsmoothcd spectrum 
to the rest- frame (2;=0), interpolating them to a common 
wavelength scale of ^O.l/zm/pix, and taking a straight 
(non-weighted) average at each wavelength position. The 
averaged spectrum is shown in Figure 2a. The most 
dominant feature is the silicate absorption feature at 

® Sio is defined in Table 3 
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TABLE 1 
Optical and IR properties. 



IRS ID MIPS Bw'^-'' R°--^ K°-^ fi/3.6'= L4.5= fi,5.8= tu^" ^24" {^70 f^lGO R-[24] 

name mag mag mag mag mjy mjy mjy mjy mjy mjy mjy mag 



BootesXl 


SST24J143048 


3- 


h322532 


24 


00 


23 


11 


22 


33 









13 


0.22 





77 


0.68 


3 


15 


<12 


<120 


14.70 


BootcsX2 


SST24J142731 


3- 


h324434 


25 


54 


23 


04 


21 


51 









09 


0.12 





20 


0.59 


2 


61 


<15 


<100 


14.42 


BootcsXS 


SST24J142609 


0- 


h333303 


23 


33 


22 


80 


21 


84 


18 


29^* 





05 


0.06 





10 


0.18 


2 


49 


<20 


<120 


14.13 


BootcsX4 


SST24J143234 


9- 


h333637 


24 


53 


22 


96 


21 


53 


18 


65^* 





18 


0.35 





73 


1.68 


2 


92 


30±6 


110±22 


14.47 


BootcsX5 


SST24J 142625 


5- 


h334051 


25 


14 


23 


53 


22 


04 


18 


OO'' 





21 


0.32 





52 


0.82 


2 


11 




<120 


14.69 


BootcsX6 


SST24J143232 


6- 


h335903 


23 


09 


22 


79 


22 


10 


17 


85^* 





04 


0.05 





16 


0.38 


2 


42 


10±2 


80±16 


14.09 


BootcsX? 


SST24J143038 


4- 


h340929 


24 


34 


22 


86 


21 


67 


18 


53^* 





06 


0.10 





17 


0.35 


2 


40 


12±2 


<100 


14.15 


BootcsXS 


SST24J143301 


5- 


h342341 


23 


48 


22 


85 


22 


54 


18 


32'' 





10 


0.16 





35 


0.62 


2 


17 


<15 


<120 


14.04 


BootcsXg 


SST24J 142849 


0- 


h343432 


23 


58 


22 


35 


21 


14 


18 


11= 





08 


0.10 





13 


0.33 


2 


36 


<15 


<100 


13.60 


BootcsXlO 


SST24,I143834 


9- 


h343839 


22 


76 


22 


44 


22 


58 









06 


-1.09 





25 


-1.09 


2 


25 


<20 


<150 


13.66 


BootcsXll 


SST24,J143715 


1- 


h345323 


24 


15 


23 


73 


23 


06 









16 


0.27 





49 


0.95 


2 


51 


<15 


<120 


15.07 


BootcsX12 


SST24,I142537 


8- 


h351735 


24 


21 


22 


12 


20 


69 


18 


30*= 





09 


0.12 





21 


0.38 


2 


28 


25±5 


<120 


13.36 


HoucklS 


SST24J143644 


2+350627 


24 


71 


23 


69 


23 


36 


19 


11/ 





04 


0.10 





32 


0.69 


2 


34 


9it3'' 


43±12'' 


14.96 



X-ray bright sub-sample 


BootesX13 


SST24 J 143723.4+334308 


24.96 


22.78 


22.49 - 0.03 


0.03 


0.04 


0.11 


1.44 


20+4 


<120 


13.52 


BootesX14 


SST24J142707.0+325213 




22.03 


0.19 


0.16 


0.19 


0.25 


1.22 


<12 


<100 


12.59 


BootesX15 


SST24J143359.0+331301 


23.95 


22.40 


21.06 17.15'' 0.23 


0.29 


0.42 


0.61 


1.92 


<15 


<120 


13.45 



^AU quoted magnitudes arc Vega magnitudes from the NDWFS DR3 (Jannuzi et al. in preparation). 
'^Errors on B-w, R, and /C-band magnitudes are <0.1. 

*^Errors on f^(3.6), fi.(4.5), f^(5.8), f^(8), and f^(24) flux densities are <0.1 mJy. 
■^Ks-mag from FLAMEX survey (Elston et al. 2006). 
^JtT-mag from NDWFS survey (Dey et al. in preparation.) 
^K'-mag from Keck/NIRC observation. 

^Because BootesXlO lies at the edge of our IRAC coverage, no detections were found in this band. 
^Value from deeper imaging data presented in Tyler et al. submitted. 



TABLE 2 



IRS ID X-ray X-ray" hardness fx*^ Iogio(fx / iopt)'^ log(Lx^) 

name counts ratio'^ (erg s~^) 



BootcsXl 


CXOXB 


,1143048 


3H 


h322531 


5+2.2 


1. 


BootesX2 


CXOXB 


J142731 


4H 


h324435 


6+2.4 


0. 


BootesXS 


CXOXB 


J142609 


OH 


h333304 


5+2.2 


0. 


BootosX4 


CXOXB 


J143234 


9H 


h333636 


4+2.0 


0. 


BootcsX5 


CXOXB 


J142625 


5H 


h334052 


7+2.6 


— 


BootcsXO 


CXOXB 


J143232 


6H 


h335902 


21+4.6 


— 


BootcsX7 


CXOXB 


J143038 


5H 


h340928 


7+2.6 


0. 


BootcsX8 


CXOXB 


J143301 


5+342343 


14+3.7 


0. 


BootcsXg 


CXOXB 


J142849 


0+343432 


6+2.4 


0. 


BootcsXlO 


CXOXB 


J143835 


OH 


h343839 


9+3.0 


0. 


BootcsXll 


CXOXB 


J143715 


IH 


h345323 


5+2.2 


— 


BootesX12 


CXOXB 


J142537 


8+351735 


29+5.4 


— 


HoucklS 


CXOXB 


J 143644 


2+350626 


7+2.6 


0. 



n lR+0.20 

Q-ia 

n 44+0.07 

72±^P6°« 
'^^0.22 
nr)+0.Ii5 

67+8:18 
0.20±J-| 

o.io+o-Sl 



8.5+5.6 


0.67 


44.22 


10.6+5.9 


0.74 


43.63 


9.5+5.9 


0.60 


44.44 


8.1+5.6 


0.59 


43.75 


14.3+6.7 


1.07 


44.23 


43.8+10.0 


1.26 


44.25 


14.0+6.7 


0.79 


43.97 


27.1+8.4 


1.07 


44.95 


12.5+6.5 


0.54 




18.0+7.3 


0.73 


44.99 


20.5+6.0 


1.30 




56.8+11.1 


1.10 




13.2+6.5 


1.10 


44.55 



X-ray bright sub-sample 

BootesX13 CXOXB J143723.5-h334307 48ib6.9 -0.21ib2-2| 92.0ibl3.9 1.58 

BootesX14 CXOXB J142707.0-h325212 34ib5.8 -0.47ibg-g| 64.9ibl2.1 1.21 

BootesXlS CXOXB J143359.1+331300 19ib4.3 0.47ibg g;^ 37.2ib9.6 1.03 



^Counts in observed total (0.5-7 keV) band and their associated Poisson errors. 

^The Hardness ratio is defined as H-S/H+S where H denotes the X-ray counts in the hard (2-7 keV) band and S denotes the X-ray counts in the 
soft (0.5-2 keV) band. 

^'Observed 0.5-7 keV X-ray flux and error in 10"^'^ ergs cm~^ s~^ assuming a power-law spectrum with photon index F— 1.7 and a galactic column 
density of 1 x 10^'^ cm^ (Stark et al. 1992). No correction has been made for intrinsic X-ray absorption. 
■^X-ray to optical (i^-band) flux ratio defined as logio(fx / fopt)— logio(fx) -\- 0.4i? -|-5.5 (Brand et al. 2006a) 
^Rest-frame 0.5-7 keV X-ray luminosity calculated using redshifts in Table 3. 
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TABLE 3 

IRS AND DERIVED PHYSICAL PROPERTIES. 



IRS ID 


e.time 


e.time 


IRS 




i/Li^8/im 




f,.(15) 


Sio'' 


PAH EW<= 


SFR9 


L/fl(SFR)'' 


IRS 




SL" (s) 


LLf (s) 


z 


log(L0) 


log(L0) 


log{Lo) 






(nm) 


(Mo yr-i) 


log(L0) 


class'* 


BootesXl 


240 


480 


1.7 


12.2 


12.6 


13.3 




-0.4 


<0.05 


<1380 


<12.8 


abs 


BootesX2 


360 


720 


0.9 


11.5 


11.8 


12.5 


3.3 


-1.0 


<0.04 


<190 


<12.0 


abs 


BootesXS 


480 


960 


2.0^" 


12.2 


12.6 


13.3 




-0.4 


<0.20 


<1170 


<12.9 


abs 


BootcsX4 


360 


720 


1.12 


12.2 


12.4 


13.1 


2.4 


-1.4 


<0.04 


<450 


<12.4 


abs 


BootcsX5 


480 


960 


1.4 


12.0 


12.3 


13.0 


2.8 


-0.4 


<0.09 


<590 


<12.5 


abs 


BootcsX6 


360 


720 


0.9 


















pow/em? 


BootesX7 


480 


960 


1.1 


11.5 


11.9 


12.6 


6.3 


-0.5 


<0.05 


<110 


<11.7 


abs 


BootesXS 


480 


960 


2.1 


12.3 


12.5 


13.2 




-0.6 


<0.05 


<560 


<12.5 


abs 


BootesXQ 


480 


960 




















pow 


BootesXlO 


480 


960 


2.6 


12.6 


12.9 


13.6 




< -1.8 


<0.10 


<1670 


<12.9 


abs 


BootesXl 1 


480 


960 




















pow 


BootesX12 


480 


960 




















pow 


Houckl3 


480 


720 


1.95 


12.3 


12.5 


13.2 




-2.9 


<0.03 


<420 


<12.3 


abs 












X-ray bright 


sub-sampl 


e 










BootesX13 


480 


1200 




















pow 


BootesX14 


480 


1200 




















pow 


BootesXlS 


480 


1200 




















pow 



^Total integration times in SLl. 

^Total integration times in each of LLl and LL2. 

'^Silicate absorption strength is defined as S^io = In f°^^^(io^) > where /obs (10/.im) is the observed flux density at the peak of the lOjum feature 
and /coiif (10/Am) is the continuum flux at the peak wavelength, extrapolated from the continuum to either side. 
*^Thc total infrared luminosity, L/i^=L(8-1000jLtm), as estimated from the measured rest-frame 8/im flux density using the relation for the template 

of Mrk 231 (L/k^S.O x /yL,48^m)) . 

*^ Absorption-dominated and featureless power-law sources are denoted "abs" and "pow" respectively, "em" denotes possible silicate emission. 
^Rest-frame equivalent width of 6.2/jm PAH emission feature. 

^The redshift could be ^^ = 1.5 if instead of silicate absorption at an observed wavelength of 29A, the spectrum is interpreted as having silicate 
emission at an observed wavelength of ^24A. However, silicate emission is generally only observed in type I AGN which tend to have much brighter 
optical magnitudes for such bright 24/jm emission (e.g.. Brown ct al. 2006). 

^Maximum star formation rate (SFR,) estimated from the measured limit on the 6.2/j-m PAH emission feature using the relation of Pope et al. 
(2007). 

^Maximum infrared luminosity associated with star formation estimated from the measured limit on the 6.2/^m PAH emission feature using the 
relation of Pope et al. (2007). 
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9.7/ini. The average silicate absorption strength (as de- 
fined in Table 3) is S'io=— 1-0 (with individual values 
ranging from —0.4 to —2.9). For comparison, we also 
show the average spectrum of the absorption-dominated 
70/xm-selected sources presented in Brand et al. (2008) 
(Figure 2b). Even though the X-ray-selected absorption- 
dominated sources have a similar average i?— [24] color to 
that of the 70/xm-selcctcd absorption-dominated sources, 
they have a much shallower average silicate absorption 
(510=— 1-0 compared to 510=— 2.1). The average sili- 
cate absorption of the X-ray sample is substantially less 
than the median value for local ULIRGs (Hao et al. 2007; 
Spoon et al. 2007) and the optically obscured sample 
observed by Houck ct al. (2005) and Wecdman et al. 
(2006b), and is slightly larger that the median value for 
Seyfert 2 AGN (-0.6; Hao et al. 2007). 

The individual spectra show no evidence for PAH emis- 
sion (6.2y^m equivalent width <0.2/im; see Table 3), al- 
though there is a hint of the 6.2 and 7.7/im PAH fea- 
tures in the average spectrum. However, because of the 
redshift uncertainties for the individual spectra, it is dif- 
ficult to be confident in any weak narrow emission line 
features that appear in the average. The possible fea- 
ture at 7.7/im could also be due to the [Ne VI] emission 
line at 7.6/im. For each individual source, we use the 
noise spectrum to estimate the maximum flux from the 
6.2/im PAH emission feature (we estimate the flux of the 
nearest noise feature to 6.2/im). We then use the rela- 
tion of Pope et al. (2007) to estimate the corresponding 
maximum SFR and infrared luminosity associated with 
star-formation. The values are shown in Table 3. In 
some cases, the spectra arc not sensitive enough to put 
a useful limit on the SFR. BootesX2 and BootesX7 have 
the lowest limits of SFR<190 Mq yr^^ and SFR<100 
M0 yr~^ respectively, corresponding to infrared lumi- 
nosities associated with star- formation of log Lir <12.0 
and <11.7 respectively. Comparing these values to the 
total infrared luminosity estimated from the rest-frame 
8/im flux density, suggests that less than 30% and 13% 
of the total infrared luminosities are from star formation 
processes. 

Weak indications are present in the average for [Ne II] , 
[Ne V], and [Ne HI] emission lines. The relative inten- 
sities of these lines are ^0.7:1:1. The presence of [Ne 
V] is a strong AGN indicator (Lutz et al. 1998; Farrah 
et al. 2007). The relative intensities show that [Ne V] is 
stronger than that of local Seyfcrts such as NGC 4151, 
Markarian 3, and CenA (Wecdman ct al. 2005). If there 
is a contribution of 12.7/im PAH emission to the [Ne II] 
line in the stacked spectrum, this would make this effect 
even larger. This may indicate a higher ionization of the 
narrow-line region by the more powerful AGN in these 
sources. 

We have measured the luminosities at 5.8/im to de- 
termine the continuum luminosity just shortward of the 
6.2/im PAH feature and at a wavelength without strong 
absorption (e.g., Brandl et al. 2006; Sajina et al. 2007; 
Hao et al. 2007: see Table 3 for individual values). The 
absorption dominated sources have median luminosity 
z/L^(6/tm) = 6 X 10*^ ergs s'^ (1.5 xlO^^ L©). These 
are very luminous sources, exceeding that of most local 
AGN in the sample of Hao et al. (2007): Sy 2 in that sam- 
ple have median z/Ly(6/im)~ 10 "Lq, and ULIRGs have 
z/L^(6/im)~ lO^^L©. Desai et al. (2007) find that local 



ULIRGs span the range i^L^(6/im) = 1x10^" - 5 x lO^^ 
L0. Our sources lie at the upper end of this range. The 
most luminous X-ray source, at i/Lj,(6/tm)=4 x 10^^ L©, 
is nearly as luminous as the most luminous source yet 
published from IRS spectra: source 9 in Houck et al. 
(2005) with i^L^(6/im)=5.9 x lO^^ L©. 

We estimate the total infrared luminosity, L//i'=L(8- 
1000/tm) of each of the absorption dominated sources. 
The conversion factors vary significantly depending on 
the template adopted for the conversion (Caputi et al. 
2007). Because we expect the infrared emission to be 
dominated by obscured AGN activity, we use the con- 
version for the template Mrk 231: L/ij=5xi/Li,(8/tm), 
where the rest-frame 8/im flux density is measured di- 
rectly from the IRS spectra. The estimates are about a 
factor of 7 lower than one would obtain using templates 
of starburst-dominated galaxies such as Arp 220, and 
this explains our smaller L/ r compared to that of Houck 
et al. (2005) for HoucklS which was estimated from 
F00183-7111. The infrared luminosities confirm that 
the absorption-dominated sources are all ultra-luminous 
(LiR >10^^ Lq) or hyper-luminous (Lir. >10^^ L©) in- 
frared galaxies. 

4.2. Power-law Sources 

Given the lack of features and/or optical follow-up 
of the power-law sources, we do not have redshifts for 
these sources. We therefore cannot construct an average 
spectrum or determine luminosities. The spectra of the 
power-law sources must have weak or absent silicate fea- 
tures (6.2/tm PAH equivalent widths ^0.2/xm) or they 
lie beyond z ~2.6, causing their absorption features to 
fall longward of the IRS wavelength limit. It is diffi- 
cult to determine whether silicate emission exists in these 
sources given our lack of redshift information and that 
the silicate emission feature is generally weaker than ab- 
sorption feature. However BootesX6 may have a weak 
silicate emission feature if it is at ^=0.9. All 3 of the 
"X-ray bright" subsample are power-law sources. 

5. DISCUSSION 
5.1. IRS spectra 

All the X-ray selected sources have infrared spectra 
consistent with their mid-IR luminosities being domi- 
nated by emission originating from AGN activity. This 
is expected given that their high X-ray luminosities im- 
ply that they host luminous AGN which should domi- 
nate their infrared spectra. We find a larger fraction of 
power-law sources (4/13 of the main sample and 7/16 
of the total (main -|- "X-ray bright") sample) than is 
found in the SWIRE sample of Weedman et al. (2006a) 
(they find only 1/9 to have a featureless power-law spec- 
trum in their sample which has an X-ray limit which is 
« 8x fainter than that of our sample). This may in part 
be due to them having accurate optical redshifts, which 
make it easier to identify features. It could also be due 
to our selection of more luminous X-ray sources, which 
are loss likely to bo heavily obscured and exhibit silicate 
absorption. This hypothesis is supported by the fact that 
6/7 of the power-law sources are also the brightest X-ray 
sources in the sample. 

To explain the lack of a silicate absorption feature in 
their IRS spectra, the power-law sources must have weak 
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Fk;. 1. — IRS spectra of X-ray selected optically faint infrared luminous galaxies ordered by their silicate absorption strengths. The IRS 
spectrum of Houckl3 was published in Houck et al. (2005). The spectra arc boxcar smoothed over a resolution element (approximately two 
pixels). In cases where we have an estimated redshift, the expected positions of the PAH emission features and silicate absorption features 
are shown. A typical noise spectrum of IRS observations with similar 24 ;um flux densities and exposure times can be found in Brand et al. 
(2008). 
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or absent silicate features or they must lie beyond z '^2.6. 
We argue that the latter is unlikely given that they are 
among the brightest X-ray sources in our sample. If at 
z >2.6, they must be extremely intrinsically luminous. 
Assuming a canonical power-law spectrum with pho- 
ton index, r=1.7, the typical X-ray flux of fo.5-7feey=4 
xlO^^^ ergs cm~^s^^ corresponds to a rest-frame X-ray 
luminosity of Lo.5_7feeV =2 xlO^^ ergss"^ at 2=2.6 which 
would suggest that they were among the most powerful 
quasars known. If at > 2.6, these sources would be sig- 
nificantly more X-ray luminous than any of the sources 
at z < 2.6. We would not expect a sharp difference in the 
luminosities of sources at higher redshifts. Furthermore, 
only 0.7% of the ~1000 X-ray sources in the Chandra 
XBootes survey (Kenter et al. 2005) with spectroscopic 
redshifts from AGES (Kochanek et al. in prep) have 
luminosities this high. In the following discussion, we 
assume that these sources are at similar redshifts to the 
rest of the X-ray selected sources and that they must 
have weak or absent silicate features. 

5.2. Infrared colors of X-ray selected AGN 

In Figure 3, we show the IRAC color-color diagram 
for our IRS sources. To increase the number of X- 
ray selected sources, we also plot the 9 X-ray selected 
IRS sources from the Lockman Hole field of the SWIRE 
survey presented by Weedman et al. (2GG6a). The 
SWIRE sample is selected to have f^(24^m)>0.9 mJy, 
fo.3-8fcey > 10~^^ ergs cm^^ s^^, and R > 22 mag. The 
X-ray flux limit still corresponds to typical AGN lumi- 
nosities at these redshifts. To determine the location in 
color-color space of these sources in comparison to the 
general 24 /im population, we also plot the distribution 
of the wlO,GGO MIPS Bootes sources with f^(24/im)>0.5 
mJy. The rough regions expected to be inhabited by 
different populations are shown. We also plot the empir- 
ical line of Stern et al. (2005) (see also Lacy et al. 2004) 
to separate AGN-dominated sources from Galactic stars 
and normal galaxies. 

20/25 of the X-ray selected sources lie within the AGN 
region of Stern et al. (2005). Hickox et al. (2007) show 
that large amounts of extinction (A^, >50) toward the ac- 
tive nucleus and/or a large contribution from a starburst 
can cause the IRAC colors of an object to fall outside of 
the Stern et al. (2005) region. Because our sources are 
largely within the AGN region, they are unlikely to be 
heavily extincted or have a large starburst contribution 
to their mid-IR emission. There appears to be no correla- 
tion between the position in IRAC color-color space with 
the presence of the silicate absorption strength, again 
supporting the hypothesis that the silicate feature arises 
in the cooler gas which is traced by longer infrared wave- 
lengths. The sources in the X-ray bright sub-sample ap- 
pear to have slightly bluer colors which may demonstrate 
that they are less dust obscured. The IRAC color-color 
classification is expected to deteriorate at z^l as the 
observed IRAC bands are pushed towards lower wave- 
lengths in which the stellar light dominates over hot dust 
emission (Stern et al. 2005; Lacy et al. 2004). However, 
it likely holds to higher redshifts for the most luminous 
sources because the near-infrared emission is still dom- 
inated by the hot dust, even in the integrated spectral 
energy distribution. 



5.3. Infrared-to- optical colors of X-ray selected AGN 

In Figure 4, we show the 24-to-8 /xm vs. 24-to-0.7 

/im color-color diagram for our sources (see Yan et al. 
(2GG4) for the location of different populations in this 
color-color space). The 24-to-8 /im and 24-to-G.7 /im 
colors are defined as f(24:8)=logio[i'fi/(24/im)/i/fi^(8/im)] 
and f(24:R)=logio[^'fi^(24/im)/!/f,^(0.7/im)] respectively. 
f(24:R) is a good indicator of obscuration and/or red- 
shift. f(24:8) has been used as a crude measure 
of the spectral slope and hence the dust tempera- 
ture distribution and whether the source is AGN- or 
starburst-dominated. Brand et al. (2GG6b) find that 
AGN-dominated sources tend to have f(24:8)<0.3 and 
starburst-dominated sources tend to have f(24:8)>0.3. 
Even though the X-ray selected sources must harbor a 
powerful AGN, their 24-to-8 /xm colors suggest that the 
sources are dominated by a combination of AGN and 
starbursts. We suggest that this diagnostic breaks down 
for obscured sources due to the heavy absorption at near- 
IR wavelengths which will steepen the infrared slope re- 
gardless of the underlying power source. In addition, at 
2 PS 1-2, the silicate absorption feature is moving through 
the 24mum filter, complicating the interpretation fur- 
ther. 

The i?-[24] colors of the X-ray selected sources are 
typically bluer than that of the most extreme IRS 
sources presented in Houck et al. (2005) and Weed- 
man et al. (2006b) which typically have R - [24] >15 
(f(24:R)>1.83). The faintest optical magnitude is rela- 
tively bright (_R=23.7) despite the fact that there is no 
optical selection criteria against including optically faint 
sources. We suggest that the X-ray-selection is biasing 
us against the most heavily absorbed (and therefore op- 
tically faintest) sources. Indeed, the sources in the X-ray 
bright sub-sample have among the bluest i?-[24] colors 
due at least partially to their relatively bright i?-band 
magnitudes. 

5.4. Correlations between Infrared and X-ray properties 

Figure 5 shows the observed X-ray flux as a function of 
X-ray hardness ratio for the combined sample of Bootes 
and SWIRE X-ray IRS sources. There are no sources 
with both small X-ray fluxes and small X-ray hardness 
ratios. The grey-scale shows that this is not seen in the 
larger population of X-ray sources in the same redshift 
range, and only becomes apparent when we jointly select 
X-ray sources that are both infrared bright and optically 
faint. Our selection criteria result in a sample of bolo- 
metrically luminous high redshift AGN. In these cases, 
small observed X-ray fluxes will be associated with more 
obscured sources which will also have larger hardness ra- 
tios. 

For our combined sample, the absorption-dominated 
sources tend to have fainter X-ray fluxes and larger X- 
ray hardness ratios (with a median value of HR=0.55) 
which are indicative of more absorbed X-ray spectra and 
larger HI column densities. Conversely, the power-law 
sources have brighter X-ray fluxes and smaller hardnc;ss 
ratios (with a median value of HR=— 0.2) indicative of 
less obscured X-ray spectra and smaller HI column den- 
sities. This suggests that the dusty absorbing regions 
responsible for the silicate absorption also contain suffi- 
cient gas to cause X-ray absorption or that both the gas 
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Fk;. 2. — Average IRS spectrum of X-ray sources with redshifts 
(top). For comparison, the average IRS spectra of absorption- 
dominated sources 70/im sources from Brand et al. (2008) is also 
shown (bottom). The expected positions of typicaUy strong emis- 
sion hues and the sihcate absorption line are shown. The average 
spectrum combines the individual un-smoothcd rest-frame spectra 
and is not smoothed. 

and dust lie along similar lines of sight. Deeper silicate 
absorption features may be associated with more em- 
bedded AGN with higher column densities of gas. Our 
results are consistent with those of Shi et al. (2006) who 
measure the silicate absorption depth for a sample of 
85 nearby AGN with X-ray data and find a correlation 
between the strength of the silicate feature and the HI 
column density. Our lack of redshifts for the power-law 
sources means that we must compare the X-ray fluxes 
and hardness ratios rather than the X-ray luminosities 
and fitted HI column densities. Any large difference in 
the redshift distribution of the two populations could af- 
fect our results (for example higher redshift sources with 
identical X-ray luminosities will have softer X-ray spec- 
tra and smaller X-ray fluxes). However, in Section 4.2, 
we argue that the power-law sources lie at z <2.6 and 
there is no reason to think that they lie at significantly 
different redshifts to the absorption-dominated sources 
given their similar selection criteria. 

Although we find a trend in the X-ray and infrared 
properties of these sources, the correlations are only 
weak. The difiiculty in establishing a clear relation be- 
tween X-ray and infrared properties may be explained 
by the distribution of their gas and dust. Risaliti et al. 
(2002) show that the X-ray absorbing column density of 



Seyfert 2 galaxies varies on yearly timescales for most 
sources. This suggests that most of the X-ray absorp- 
tion must occur on sub-parsec scales. They obtain even 
tighter constraints in Risaliti et al. (2007) for a Seyfert 
galaxy which exhibits spectral variation on timescales 
of days. This suggests that, at least for less powerful 
sources, the obscming gas is on a similar scale to that 
of the broad-line region. Significant amounts of dust are 
unlikely to lie at these distances because of the dust sub- 
limation temperatures. The difference in the location of 
much of the gas and dust may explain the large scat- 
ter between parameters tracing the X-ray and mid-IR 
obscuration (e.g., Rigby et al. 2004). Lcvenson et al. 
(2007) suggest that the obscuring material in dusty AGN 
is likely to consist of both dense clumpy clouds of dust 
and a smoother inter-cloud medium. They argue that 
deep silicate absorption can only be produced in a source 
embedded in a smooth medium since clumpy clouds will 
provide direct views of portions of illuminated clouds 
which will fill in the resultant silicate absorption feature. 
Since our X-ray selected sources have shallower silicate 
absorption features than typical ULIRGs, this may sug- 
gest that they are dominated by clumpy obscuration. Al- 
ternatively, they may simply have a smaller line of sight 
obscuration. 

Sturm et al. (2006) obtained IRS spectra for a sample 
of 7 X-ray selected AGN which were specifically chosen to 
have large X-ray luminosities and large column densities 
of gas (due to their large X-ray hardness ratios). Given 
the correlation between HI column density and silicate 
absorption strength (see discussion above and Shi et al. 
2006), we might expect them to exhibit deep silicate fea- 
tures. However, all of their sources exhibit featureless 
power-law spectra. The solution to this apparent prob- 
lem may be that they selected a highly unusual popula- 
tion of sources which were both extremely X-ray lumi- 
nous and X-ray obscured. Figure 5 shows that although 
X-ray obscured sources are more likely to have silicate 
absorption features, X-ray luminous sources are likely to 
have featureless power-law spectra. Perhaps more lu- 
minous AGN are more capable of breaking up the sur- 
rounding dust into discrete clouds which would fill in the 
silicate absorption feature in these cases despite the large 
column densities of gas. 

We compare the fraction of power-law dominated 
to absorption dominated sources in different samples. 
Weedman et al. (2006a) probe the lowest X-ray lumi- 
nosities and find only i/9 (11%) are power-law domi- 
nated. Our sample has higher X-ray luminosities and we 
find the fraction of power-law to absorption dominated 
sources to be 4/13 (31%) in the main sample and 7/16 
(43%) in the total (main + X-ray bright) sample. All of 
the sources in the Sturm ct al. (2006) study are power- 
law dominated. Their sources are particularly bright in 
the X-ray, with intrinsic X-ray luminosities in the range 
logLa;=44.0-45.5. This suggests that more X-ray lumi- 
nous sources are more likely to have shallower silicate 
absorption (which will appear as power-law sources). 

Krabbe et al. (2001), Lutz et al. (2004), and Horst 
et al. (2006) show that there is a strong relation be- 
tween AGN hard X-ray emission and mid-IR continuum 
in low redshift active galaxies. Figure 6 shows the AGN 
continuum luminosity at 6/Km as a function of 2-7 keV 
X-ray luminosity for the absorption-dominated sources 
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Fig. 3. — IRAC color-color diagram for X-ray selected IRS sources (large symbols). The Bootes X-ray selected sources are denoted 
by large blue symbols (with the largest symbols denoting the X-ray bright sub-sample). The fainter X-ray sources from Weedman et al. 
(2006a) are also shown (smaller red symbols). The sources are split into silicate absorption-dominated sources (downwards triangles), and 
power-law dominated sources (stars). We show the distribution of all ailO,000 MIPS sources with f^(24/im)>0.5 mjy in greyscale. The 
rough regions expected to be inhabited by different populations are also shown. The black dashed line is that proposed by Stern et al. 
(2005) to empirically separate AGN from Galactic stars and normal galaxies. The arrow represents the reddening vector, estimated using 
the extinction curve of Draine (2003), for a source at z=l reddened by a dust column corresponding to Ru=3.1. 
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Fig. 4. — 24-to-8 fim vs. 24-to-0.7 fim color-color diagram for X-ray selected IRS sources. The 24-to-8 /im and 24-to-0.7 /im colors 
are defined as f(24:8)=logio[i^fiy(24^m)/!^f^(8^m)] and f(24:R)=logio[i^fiy(24^m)/!^fy(0.7/im)] respectively. The symbols and colors are the 
same as in Figure 3. The vertical dashed line shows the R-[24]>14 (log(!^f^ (24)/i^f^ (R)>1.43) criteria used to select the powerful obscured 
sources presented in Houck et al. (2005) and Weedman et al. (2006b). The horizontal dashed line shows the 24-to-8/im color criteria used 
by Brand et al. (2006b) to roughly divide steeper spectrum starburst sources from shallow spectrum AGN-dominated sources at z >0.6. 



Spitzer IRS Spectra of X-ray Sources 



11 



-0.5 - 
-1.oLl 




-15.5 -15.0 -14.5 -14.0 -13.5 -13.0 
Log (0.5-7 keV X-ray Flux) (ergs cm"^ s ') 

Fig. 5. — 0.5-7 keV band X-ray flux against X-ray hardness ratio 
for X-ray selected IRS sources. The hardness ratios can range from 
1 (hardest X-ray spectrum) to -1 (softest X-ray spectrum). The 
Bootes X-ray selected sources are denoted by large blue symbols 
(with the largest symbols denoting the X-ray bright sub-sample). 
The fainter X-ray sources from Weedman et al. (2006a) are also 
shown (smaller red symbols). The sources are split into silicate 
absorption-dominated sources (downwards triangles), and power- 
law dominated sources (stars). Poisson error bars are also shown. 
The grey scale shows the distribution of all X-ray sources in the 
XBootes survey with redshifts z=0.9-2.6. 

(the power-law sources could not be included because we 
do not have redshifts to derive the luminosities). Al- 
though the absorption-dominated sources do not show a 
clear correlation between these two quantities, this lack 
is likely to be due to the large luminosity errors and rel- 
atively narrow luminosity range of the sample. When 
we compare our sources with lower luminosity sources at 
low redshift (Lutz et al. 2004), we find that they form 
the higher luminosity extension of the relation. We have 
chosen not to correct the X-ray luminosity for absorp- 
tion given the large uncertainties. 7/9 of our sources lie 
at L(2-7keV)/i/L^(6^ni)<0.1, compared to only 8/47 of 
the Lutz et al. (2004) sources (before correcting for ab- 
sorption). This suggests that our sources are generally 
more obscured than the low redshift, lower luminosity 
AGN presented in Lutz et al. (2004). 



6. CONCLUSIONS 

We have presented mid-infrared spectra obtained us- 
ing Spitzer IRS of a sample of 16 X-ray loud 24 fim 
sources. The spectra of these sources are character- 
ized by either silicate absorption or a featureless power- 
law, with no evidence of PAH emission. This suggests 
that the mid-infrared emission from these sources is 
dominated by warm dust emission surrounding the ac- 
tive nucleus with little contribution from star forma- 
tion. The absorption-dominated sources lie in the range 
0.9< z <2.6 and have very high X-ray and infrared lumi- 
nosities (log L/i^=12.5-13.6 Lq). The average absorption 
depth in our absorption-dominated X-ray sources is simi- 
lar to that of local Seyfert 2 AGN and much smaller than 
that of typical ULIRGs. This suggests a difference in the 
distribution of the dusty material between our sample 
and typical ULIRGs. The weak silicate absorption fea- 
ture in our sample could indicate a more clumpy dust 
distribution or a lower line of sight obscuration. There 
is evidence in the average spectrum for weak [Ne II] , [Ne 
V], and [Nelll] lines. The [Ne V] is relatively stronger 
than is seen in local Seyferts, suggesting a higher ion- 
ization of the narrow-line region in these more powerful 
AGN. Given their higher X-ray fluxes, it is unlikely that 
the power-law sources lie at z >2.6. Instead, they must 
have only weak silicate absorption or emission features. 
Comparison with other samples suggests that the frac- 
tion of power-law sources appear to increase with X-ray 
luminosity and decrease with X-ray obscuration. Both 
the infrared and infrared to optical colors suggest that 
the sources are relatively unobscured compared with R- 
[24]>15 sources. Our X-ray selection is likely to be pref- 
erentially biased to less obscured and/or more powerful 
AGN. We find that sources with silicate absorption fea- 
tures tend to have fainter X-ray fluxes and larger hard- 
ness ratios. This suggests that more absorbed sources 
also have higher column densities of absorbing gas along 
the same line of sight. 




41 42 43 44 45 46 
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Fig. 6. — Rest-frame 6fim luminosity against rest-frame 
absorption-uncorrected 2-7 KeV X-ray luminosity for absorption- 
dominated Bootes sources (blue triangles). The Poisson error bars 
are shown for the X-ray luminosity. The low redshift ISO sources 
from Lutz et al. (2004) are plotted as small squares (connected 
empty and filled squares denote X-ray luminosities uncorrected 
and corrected for X-ray absorption respectively). The solid line 
shows a direct correspondence between the X-ray and Gfira lumi- 
nosities. The dotted lines show where L(2-7keV)=0.1 t'L^(6/im) 
and L(2-7keV)=0.01 uL^(6fim). 
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